Imaginary-time method for radiative capture reaction rate 
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We propose a new computational method for astrophysical reaction rate of radiative capture 
process. In the method, an evolution of a wave function is calculated along the imaginary-time axis 
which is identified as the inverse temperature. It enables direct evaluation of reaction rate as a 
function of temperature without solving any scattering problem. The method is tested for two-body 
radiative capture reaction, 16 0(a, 7) 20 Ne, showing that it gives identical results to that calculated 
by the ordinary procedure. The new method will be suited for calculation of triple-alpha radiative 
capture rate for which an explicit construction of the scattering solution is difficult. 



I. INTRODUCTION 

Radiative capture reaction rate far below the Coulomb 
barrier is an essential input for quantitative understand- 
ing of stellar evolution and nucleosynthesis However, 
direct experimental measurements of relevant cross sec- 
tions far below the Coulomb barrier often accompany 
difficulties because of their exponentially small cross sec- 
tions. There are also a few three-body processes of sig- 
nificance for which experimental measurements are not 
feasible. Theoretical evaluation of the radiative capture 
reaction rate is thus important. 

The radiative capture rate is composed of two distinct 
contributions, resonant and nonresonant processes. A 
potential model is often employed for the theoretical eval- 
uation of the nonresonant capture rate. For two-body ra- 
diative capture processes, it is a routine procedure once 
the model potential is given. One first solves the radial 
Schrodinger equation (coupled channel equation if nec- 
essary) under an appropriate scattering boundary condi- 
tion to obtain scattering cross section. One then calcu- 
lates the capture reaction rate as a function of temper- 
ature by integrating the cross section over the collision 
energy with an appropriate Boltzmann weight. 

Theoretical evaluation of radiative capture reaction 
rate for three-body processes is a much more difficult 
problem. It is well recognized that the triple-alpha ra- 
diative capture reaction to form 12 C is a key process to 
produce heavy elements 0. At a temperature above 0.1 
GK, a resonance state of 0^ of 12 C which is known as 
the Hoyle state contributes dominantly. Below 0.1 

GK, on the other hand, nonresonant contribution is con- 
sidered to be significant 0-0] • Recently, Ogata ct.al. Q 
conducted a serious evaluation of the rate with the CDCC 
(Continuums-Discretized Coupled-Channels) method, a 
three-body reaction theory which has been successful for 
nuclear direct reactions [3, ITo| . The radiative capture 
rate which Ogata et al. reported was surprisingly large 
below 0.1 GK in comparison with the rate which has been 
employed in standard steller evolution calculations [l[. 

Theoretical evaluation of three-body radiative capture 
rate accompanies several difficulties. It is by no means 



obvious how to define theoretically the cross section of 
the triple-alpha radiative capture process, because an an- 
alytic asymptotic form of the scattering wave function of 
three chargcd-particles is not known. One also needs to 
solve the three-body problem in a huge spatial region for 
reactions far below the Coulomb barrier, since the alpha 
particles pass through a barrier for a long distance to 
penetrate it. 

In Ref. [ll|, de Diego ct.al. proposed an alternative 
procedure for the calculation of triple-alpha capture rate. 
They consider an inverse process, a photo-absorption of 
12 C in the excited 2 + state, and calculate the transition 
probability in the bound state approximation. This pro- 
cedure allows one to avoid the difficulty of calculating 
scattering solution for three chargcd-particles. However, 
a number of bound states need to be calculated in their 
approach. 

In this paper, we propose a new computational method 
for the radiative capture rate. We will show that the ra- 
diative capture rate as a function of temperature may be 
calculated directly by solving an equation which looks 
like a time-dependent Schrodinger equation along the 
imaginary-time axis. The new method requires neither 
any solutions of scattering problem nor any bound state 
solutions except for a final bound state wave function 
after the capture. Since the new method allows us to 
avoid the difficulties mentioned above, we consider it will 
be useful for the calculation of radiative capture rate of 
triple-alpha process. In this paper, we demonstrate fea- 
sibility of the method by applying it to two-body capture 
reaction, 16 0(a, 7) 20 Ne, as an example. It will be shown 
that the new method gives an identical result to that 
calculated by the ordinary method using the two-body 
scattering solution. 

The construction of this paper is as follows. In 
Sec. Ill Al we present the imaginary-time formalism for 
the radiative capture reaction rate. In Sec. Ill B[ we dis- 
cuss how the resonant and nonresonant contributions are 
included in our formalism. In Sec. IIIII our method is 
exemplified by applying it to the two-body capture reac- 
tion, 16 O(a,7) 20 Nc. In Sec. IIII Al we summarize radial 
equations to be employed in the practical calculation. In 
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Sec. 1111 B[ we show results with ordinary method solving 
two-body scattering problem. We then compare the re- 
sult of the new method with that of ordinary method in 
Sec. IIII Cl Sec. [TV] is devoted to summary. 



II. THEORY 

A. Imaginary-time method for radiative capture 
rate 

We consider a radiative capture process of two or three 
nuclei confined in a large spatial area of volume V. The 
transition rate between nuclear states i and / accompa- 
nying an emission of a photon of multipolarity Xfj, is given 
by H3, 
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where is a transition operator. The energy of emit- 
ted photon E~f is equal to the energy difference of two 
states, E 1 = Ei — Ef. The initial state ^ is a scatter- 
ing state in which i specifies the relative momentum of 
colliding nuclei and other quantum numbers. The final 
state i$> t is a bound state wave function after emitting 
the photon. 

For two-body collisions, transition rate in a unit spa- 
tial area and in a unit time under unit number densities 
of colliding nuclei is given by VTj^ and is equal to ucr/i, 
where v is the relative velocity and ay, is the cross sec- 
tion. For three-body collisions, the reaction rate in a unit 
spatial area and in a unit time is given by V 2 T^^\ 

We denote the inverse temperature as j3 = l/fc^T and 
express the thermonuclear reaction rate at the inverse 
temperature j3 as r{0). This is related to the transition 
rate by 
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where N = 2 for two-body and N = 3 for three-body col- 
lisions, respectively. Mf indicates the magnetic quantum 
number of final state /. The denominator is evaluated to 
be 
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where \i is the reduced mass, [i = ra\m-2.livri\ + m<i) for 
two-body case and fi = mxrnafns f {rrii + m,2 + m^) for 
three-body case. u>i accounts for the degeneracy of the 
initial state. 

An essential trick which brings us an imaginary-time 
evolution formula for the reaction rate is an employment 
of the spectral representation of the Hamiltonian. Let 



f(H) be a certain function of the Hamiltonian operator 
H . We then have 

f(H) = /(£»)|*»><*n|+ E f(Ei)\*i)(*i\ 

unbound i£ scattering 

(4) 

where E n and 4>„ are energy eigenvalues and eigenfunc- 
tions of bound states and Ei and $j are those of scatter- 
ing states. For a two-body scattering state, the energy Ei 
is given by Ei = h 2 k 2 /2^, where k specifies the relative 
wave number of colliding nuclei. 



Employing Eq. ((4]) with f(x) = e 
may rewrite Eq. ^ as 
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where P is a projector to remove bound states, 

P=l- ]T l*n)(*n|- (6) 

unbound 

Equation (J5]) is the principal result of this paper. We 
find that the initial scattering states are removed in this 
expression. 

For a practical calculation of Eq. (|5|), we introduce a 
wave function by 
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Then the reaction rate is expressed as 

8tt(A + 1) 



1 / 2^/3?^ 3/2 



to, V ft J ^A((2A + 1)!!) 2 
E(*/|M Am |* Am ,/(/3)}. 



(8) 
(9) 



M fli 



The wave function ^\^j((3) satisfies a time-dependent 
Schrodinger equation along the imaginary-time axis, 



with the initial condition, 
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In practical calculations, we repeat evolutions with a 
small imaginary-time step A/3 to achieve a finite evolu- 
tion, 



(nA/3) = Pe- A ^P ■ ■ ■ Pe~ Af3k ^ Ml j(0). (12) 
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The operation of the evolution operator with a small 
imaginary-time step, e~ A/3H may be achieved with the 
Taylor expansion method, 



0> w (/3 + A/3) = Pe- A ^ AM>/ (/3) 
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In an analytic expression, the projector P is necessary 
only once in Eq. ([5]), since the Hamiltonian H commutes 
with the projector P. In practical calculations, however, 
it is indispensable to apply the projector at each step of 
Eq. HU). 



where Mr specifics a magnetic substate of the resonance. 
Noting that the perturbation theory gives an expression 
for the radiative decay width of the resonant state $r as 
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we arrive at the following result. 
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B. Resonant and nonresonant contributions 

In the ordinary treatment of radiative capture pro- 
cesses, contributions of sharp resonances are treated sep- 
arately from the nonresonant contribution. For a two- 
body collision, a contribution of the resonance of energy 
Er and width V to the reaction rate is given by [l[ 
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where T,; and Tf are partial widths of the resonance to 
the initial channel through barrier penetration and to the 
final state through 7 emission, lor is the statistical factor 
given by 

W *-(2/ 1 + l)(2/ 2 + l)' (15) 

where Jr is the spin of the resonance and A(2) is the spin 
of colliding nucleus 1(2). 

Equation ([5]) includes both resonant and nonresonant 
contributions since all the final states are summed up. To 
confirm that resonant contribution is included in Eq. ([3]) , 
we show below that the resonant contribution r^(/3) may 
be extracted from it. 

We assume that the partial decay width for gamma 
emission, T^, is much smaller than the partial width for 
binary or ternary decay through Coulomb barrier, IV 
Indeed, this is the condition that we may start with the 
transition rate expression of Eq. ^ in perturbation the- 
ory. We thus assume that the partial width decaying into 
the initial channel, I\, almost exhausts the total width, 
r ~ Tj. For a sharp resonance, we may express the reso- 
nant state by a normalized wave function &r. To calcu- 
late the resonant contribution, we replace the projector 
P in Eq. §5§ with the projector of the resonant state, 
^fl)(^ > i?|- Then we find the contribution of the resonant 
state may be expressed as 

3/2 



1 



/' 



8tt(A + 1) 
ftA((2A + l)!!) 2 ' 



-/3E R 



' Er — Ef 

he 



2A+1 



M tl j,M R 



} \M x ^r)\\ (16) 



This is equal to r^(/3) if we assume T = T, in Eq. (|T4|) . 

In practical calculations, there are two options when a 
sharp resonance exists. One is to treat the resonant con- 
tribution separately employing Eq. (|14[) . removing the 
resonant contribution from the expression of Eq. (0 by 
adding the projector of the resonant state to the projec- 
tor P. The other is just to perform the imaginary-time 
calculation as it is, so that the resonant contribution is 
automatically included in Eq. (| 14[) . 



III. TEST CALCULATION : 16 O(a,7) 20 NE 
CAPTURE RATE 

A. A potential model and radial equations 

To confirm that the imaginary-time method explained 
in the previous section works in practice, we apply the 
method to a radiative capture process of two-body colli- 
sion, 16 0(a, 7) 20 Ne. We assume a simple potential model 
for the initial a- 16 scattering state and for the final 
excited state of 20 Ne. This potential model has been 
adopted in [l~3| and has been shown to describe the pro- 
cess reasonably. 

Numerical calculations will be achieved in the partial 
wave expansion. We first summarize the formula in the 
partial wave expansion for a- 16 collision. We introduce 
a radial wave function for the relative motion in the or- 
dinary way. For bound states, we denote 
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where n is the nodal quantum number. We assume a 
normalization relation J dr\u n i(r)\ 2 = 1 as usual. For 
scattering states, we denote the radial wave function of 
energy E as uei{t) for which we assume the following 
normalization in the asymptotic region, 



with E = h 2 k 2 /2fi. Then, there follows the following 
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completeness relation for each I value, 



^] u n i(r)u nl (r') + / dEu El (r)u E i{r') = S(r - r'). 



(22) 

In the ordinary method, we first calculate the cross 
section and then calculate the reaction rate by integrating 
the cross section with a Boltzmann weight and a photon 
phase space factor. Denoting the radial wave function of 
initial state by u^. (r) and that of final state by (r), 
the reaction rate is given by 
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where we introduced QijXP) by 



dEe 



-0E 



E-E f 



2A+1 



dru 



(/) 



(24) 



In the imaginary-time method, we employ a spectral 
representation of the Hamiltonian to remove the scatter- 
ing wave function (r) . The final expression written 
in terms of the radial wave function is given by 
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where Hi i is the radial Hamiltonian for the partial wave 
li. The Pi i is the radial projector to remove bound states 
of the partial wave U. 

Introducing a radial wave function u^J. (r, (3) by 
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Namely, we consider below the case of A = 2, — 0, and 
l f = 2. 

We assume a simple Woods-Saxon form for the a- 16 
potential with a radius parameter i?o = 2.72 fm and a 
diffuseness parameter a = 0.85 fm. The depth of the 
potential is so chosen that the energies of bound states 
are reproduced reasonably. The potential depth of I = 
channel is set to Vq = —150.23 MeV to reproduce the 
ground state energy of 20 Ne from the a- 16 threshold, 
—4.63 MeV. The potential depth of / = 2 channel is set 
to Vq = —147.95 MeV to reproduce the excitation energy 
of first 2+ state of 20 Ne, 1.63 MeV. There appear many 
bound states in this potential besides the above physical 
states. They correspond to the Pauli-forbidden states 
of many-body wave function. We include all the bound 
states, both physical and Pauli-forbidden states, in the 
projection operator Pi i . The calculations shown below 
are achieved with a radial grid of Ar = 0.1 fm. The 
Runge-Kutta method is used to solve the radial equation 
from the origin and a simple five-point finite-difference 
formula is used in the imaginary-time evolution for the 
second-order derivative operator in . 



B. Ordinary method 

Before showing results with the imaginary-time 
method, we first show calculations in the ordinary ap- 
proach solving the radial Schrodinger equation for each 
incident energy. Figure [IJa) shows the radial wave func- 
tions of initial and final states. The initial scattering 
wave, u^). = o( r )' is shown by dashed curve. The in- 
cident relative energy is set to E = 0.1 MeV, which 
approximately corresponds to the Gamow window en- 
, ergy at T = 10 7 K. The final bound-state wave function, 
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The radiative capture reaction of 16 O(a,7) 20 Ne pro- 
ceeds dominantly from the initial a- 16 scattering state 
with li = relative angular momentum to the final 2 + 
state of 20 Ne at the excitation energy of 1.63 MeV after 
emission of E2 gamma ray [l3j . Since our purpose here 
is to show the applicability of our new method, we con- 
centrate on the calculation of this transition component. 
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=2 (r), is shown by solid curve. 



Figure [TJb) shows the integrand of the radial matrix 

element appearing in Eq. (j24|) . u~^i s= 2{ r ) r2u ^Eii= o( r )- ^ s 
seen from the figure, a dominant contribution comes from 
a spatial region where the final wave function u£ \ =2v) 
decays exponentially. We find the radial integration up 
to 30 fm in Eq. (|24|) is required to obtain a fully converged 
result. 

We show in Fig. [5] the reaction rate q\ ~ 2 l _o(P) °f 
Eq. (|24[) as a function of inverse temperature /3 = 1/ksT 
by solid line (left scale) . We also show the Gamow energy 
as a function of inverse temperature by dashed line (right 
scale). The inverse temperature /? = 10 MeV -1 corre- 
sponds approximately toT=l.lxl0 9 K and /3 = 1000 
MeV" 1 to T = 1.1 x 10 7 K. 



C. Imaginary-time method 
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We show in Fig. [3] the reaction rates, q\ \ ~ 2 ; j= oj calcu- 
lated with the imaginary-time method. We find the cal- 
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FIG. 1: (a): The radial wave function of the initial scatter- 
ing state, „(r), is shown for the incident relative energy 
E — 0.1 MeV by dotted curve, and the radial wave func- 
tion of the final state, _ 2 (t'), is shown by solid curve. 

(b): The overlap function, ui'^^M 7 " 2 '"^; =o( r )i appearing 
in the integrand of Eq. (|24|1 in the text. 
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FIG. 2: Solid curve and left scale: radiative capture reaction 
rate, ^1^2^1=0(^)1 calculated with the ordinary method of 
Eq. (|24p. Dashed curve and right scale: Energy of Gamow 
window as a function of inverse temperature /3. 
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culated reaction rate depends on the radial region where 
the imaginary-time evolution is calculated. In Fig. [3l 
reaction rates calculated with different choices of radial 
cutoff distance R m ax are compared. The reaction rate in 
the ordinary method, which was shown in Fig. [21 is also 
shown for comparison. 

As is seen in the figure, the reaction rate falls off too 
rapidly if the radial cutoff distance is not sufficiently 
large. The reaction rate calculated with the radial cut- 
off distance R rna x = 500 fm coincides almost completely 
with the reaction rate calculated in the ordinary method 
for a whole temperature region shown in the figure. We 
thus conclude that, to obtain accurate reaction rate at 
low temperature with the imaginary-time method, one 
needs to calculate the imaginary-time evolution of the 
wave function in a sufficiently large radial space, up to 
500 fm for T ~ 10 7 K. 

Figure [3] indicates that the reaction rates by the 



FIG. 3: The radiative capture reaction rate '_%\ . =0 (/3) cal- 
culated by the imaginary-time method of Eq. (|25[) is shown for 
several choices of radial cutoff distance Rmax ■ For a compari- 
son, the reaction rate calculated with the ordinary method of 
Eq. Q24p is also shown, which is denoted as "ordinary". 



imaginary-time method decrease exponentially at large 
f3 when the radial cutoff distance is not sufficiently large. 
We express the asymptotic behavior as qjZ% i.—qW) — 
e~P e , where the slope constant e depends on the radial 
cutoff distance, R m ax- The e increases as the radial cut- 
off distance decreases. In the imaginary-time calcula- 
tion, the wave function u^J, (r, /3) is dominated by the 
eigenfunction of the lowest eigenvalue when ft is suffi- 
ciently large. Since all the bound states are removed in 
the imaginary-time evolution by the projection operator, 
the slope parameter e should coincide with the lowest 
positive energy-eigenvalue of the Hamiltonian. Since the 
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FIG. 4: The a- le O potential is denoted by a red solid curve 
(left scale). Five lowest positive-energy eigenfunctions of 
the potential are shown by dotted curves. Calculations are 
achieved in the radial region up to 500 fm. 



Coulomb potential decreases monotonically as a function 
of radial coordinate, the eigenfunction of the lowest pos- 
itive eigenvalue should localize in the region close to the 
radial cutoff distance, if there is not a sharp resonant 
state below that energy. 

In Fig. 01 we show the a- 16 potential and the eigen- 
functions belonging to several positive low-energy eigen- 
values. Calculation is achieved in the radial region up to 
Rmax = 500 fm. The a- 16 potential is composed of nu- 
clear (Vn) and Coulomb (Vc) potentials, and the lowest 
positive-energy eigenvalue is close to the minimum of the 
Coulomb potential energy at the radial cutoff distance, 
E min ~ Z 1 Z 2 e 2 /R max . For R max = 500 fm, the energy 
is E mm ~ 0.046 MeV. For R max = 100 and 200 fm, 
E m in ~ 0.23 and 0.115 MeV, respectively. These values 
explain the slope of the reaction rate at large j3 seen in 

Fig. ED 

It is evident that the imaginary-time evolution in the 
radial region inside a certain radial cutoff distance Rmax 
takes only into account the tunneling process of energy 
higher than e 2 Z\Z^J R max - As the temperature becomes 
lower, one needs to calculate the reaction rates in a wider 
radial region. We may estimate the radial cutoff distance 
Rmax to obtain a reliable reaction rate for a given tem- 
perature /3 considering the energy of the Gamow window. 
Employing a standard formula for the peak energy of the 
Gamow window as a function of temperature and equat- 
ing the energy with the Coulomb potential energy at the 
radial cutoff distance, we obtain 
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This gives R max = 85 fm for f3 = 10 2 MeV" 1 and R max = 
394 fm for (3 = 10 3 MeV" 1 . This estimation coincides 
with the observation in Fig. [3] that the calculation up 
to Rmax = 100 fm describes reaction rate for j3 < 100 
MeV" 1 and the calculation up to Rmax = 500 fm for 
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FIG. 5: Imaginary-time evolution of the wave function as 
a function of radial coordinate r, itj^Z^Y^o^"^) defined by 
Eq. (T26]), ( a ) at /3 = 0, 0.001, 0.01, and 1 MeV" 1 in the linear 
scale and (b) at P = 100, 500, 1000 MeV" 1 in the logarithmic 
scale. 



j3 < 1000 fm. 

For a deeper understanding of the imaginary-time 
method, we show in Fig. EJa) (b) the evolution of the ra- 
dial wave function u[ =2 i 4 =o( r >^) m ^ nc imaginary-time 
for several values of ft. The calculation is achieved with 
Rmax = 500 fm. In the top four panels, the wave func- 
tions of j8 < 1 MeV" 1 are shown in the linear scale. In 
the bottom panel, the absolute value of the wave func- 
tions are shown for large /3 values in logarithmic scale. 

To start the calculation, we prepare the radial wave 
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function of the final state u„ \ _ 2 (r) inside a region, 

Rmax = 30 fm. In the top panel of Fig. EJa), the wave 
function at j3 = is shown. We find a number of spikes 
in the initial wave function ^[^"^.^(r, j3 = 0). In par- 
ticular, an intense spike is seen at around r ~ Rm!x = 30 
fm. These sharp structures originate from the operation 
of (H^ ~Ef) 5 in preparing the wave function at (3 = 0. It 
works to emphasize components with high energy eigen- 
values of the radial Hamiltonian. 

At first sight, the existence of such sharp struc- 
tures looks unfavorable and problematic. However, 
these spikes disappear immediately after we start the 
imaginary-time evolution and they will not affect the 
reaction rate at low temperature. Even at (3 = 0.001 
MeV -1 , these spikes are substantially reduced. They dis- 
appear almost completely at f3 = 0.01. 

As the inverse temperature /3 increases, the wave func- 
tion starts to shift outwards. At f3 = 1 MeV" 1 , the 
amplitude of the wave function shows a peak at around 
6 fm. As seen in Fig. [5jb), the dominant component 
of the wave function gradually shifts towards a region of 
large radial distance. Eventually, at f3 > 500 MeV" 1 , the 
radial wave function is dominated in the region of large 
radial distance. At /3 = 1000 MeV" 1 , the wave function 
of small radial region (r < 30 fm), which contributes to 
the radiative capture rate, is much smaller than that in 
the asymptotic region by about 10 20 order of magnitude. 
Thus, the imaginary-time calculation should be achieved 
with high accuracy to describe 10 20 difference of magni- 
tude of the wave function in different radial regions. 

To confirm that the result does not depend on the ra- 
dial region in which we prepare the final wave function, 
we compare reaction rates employing final wave functions 
u nf i f ( r ) prepared in the radial region with different cut- 
off radius, Rmax- 

In Fig. [5J we compare reaction rates calculated by the 
imaginary-time method employing final wave functions 
of different radial cutoff distance, Rmax- We find the 
calculated reaction rate is quite insensitive to the radial 
cutoff distance. We thus confirm that a number of sharp 
structures seen in the top panel of Fig. (Sfa), especially 

prominent at around the radial cutoff distance, Rrllx: do 
not have any influence on the reaction rate calculation. 
As seen in the inset, the reaction rate is convergent if 
we choose Rmax > 25 fm, which is consistent with our 
observation in Sec. IIII Al 

We finally present a note on the dependence of the 
reaction rate on the choice of the nuclear potential in 
the imaginary-time evolution. In Fig. [71 we compare the 
reaction rates changing the depths of the nuclear poten- 
tial, Vo of the Woods-Saxon potential in the initial scat- 
tering state with li = 0. All the other parameters are 
set to be the same. As seen from the figure, the reac- 
tion rate is quite insensitive to the choice of the param- 
eter Vq- Even without the nuclear potential, i.e. with 
Vo = 0, the reaction rate is given almost correctly. As 




icr 100 ' 1 1 1 1 1 1 1 1 1 1 

100 200 300 400 500 600 700 800 900 1000 
j3[1/MeV] 

FIG. 6: The dependence of the reaction rate gj ; =2 V=o iP) on 
the radial cutoff distance, Rill x , in preparing the final-state 
wave function v-l?Ji f= 2( r )- 




IQ- 100 I 1 1 1 1 1 1 1 1 1 1 

100 200 300 400 500 600 700 800 900 1000 
j3[1/MeV] 

FIG. 7: Reaction rates with different nuclear potential in 
the initial scattering channel are shown. The depths of the 
Woods-Saxon potential, Vo is varied. 



the inset shows, the difference is within a factor of 1.5 
in 950 MeV" 1 < /? < 1000 MeV" 1 . We thus conclude 
that the nuclear potential in the initial channel, which 
will be used in the imaginary-time evolution, has very 
small effect on the reaction rate. Of course, this conclu- 
sion applies only to the nonresonant contribution. The 
resonance energy and width are sensitive to the nuclear 
potential, and so is the resonant contribution to the re- 
action rate. 
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IV. SUMMARY 

In this paper, we proposed a new computational 
method for radiative capture reaction rate. Employing 
a spectral representation of the Hamiltonian, we have 
shown that the reaction rate as a function of temper- 
ature may be calculated without solving any scattering 
problem. Starting with an initial wave function which 
includes the final bound-state wave function after the 
emission of photon, the reaction rate as a function of 
inverse temperature, /?(= 1/fcsT), can be obtained di- 
rectly by solving a time-dependent Schrodinger equation 



in the imaginary-time axis. 

To show feasibility of the method, we show applica- 
tion of the method to 16 O(a,7) 20 Ne reaction in a sim- 
ple potential model. We have confirmed that the new 
method gives an accurate reaction rate if we solve the 
imaginary-time evolution equation in a sufficiently large 
spatial area. Since the new method does not require any 
solution of scattering equation, it will be a promising 
approach for the reaction rate of triple-alpha radiative 
capture process. The application to that process is now 
in progress. 
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